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ABSTRACT 

We analyse deep images from the VISTA survey of the Magellanic Clouds in the 
YJKs filters, covering 14 deg^ (10 tiles), split into 120 subregions, and compris¬ 
ing the main body and Wing of the Small Magellanic Cloud (SMC). We apply a 
colour-magnitude diagram reconstruction method that returns their best-fitting star 
formation rate SFR(t), age-metallicity relation (AMR), distance and mean reddening, 
together with 68% confidence intervals. The distance data can be approximated by a 
plane tilted in the East-West direction with a mean inclination of 39°, although devi¬ 
ations of up to ±3 kpc suggest a distorted and warped disk. After assigning to every 
observed star a probability of belonging to a given age-metallicity interval, we build 
high-resolution population maps. These dramatically reveal the flocculent nature of 
the young star-forming regions and the nearly smooth features traced by older stellar 
generations. They document the formation of the SMC Wing at ages <0.2 Gyr and 
the peak of star formation in the SMC Bar at ^40 Myr. We clearly detect periods of 
enhanced star formation at 1.5 Gyr and 5 Gyr. The former is possibly related to a new 
feature found in the AMR, which suggests ingestion of metal-poor gas at ages slightly 
larger than 1 Gyr. The latter constitutes a major period of stellar mass formation. We 
conhrm that the SFR(t) was moderately low at even older ages. 

Key words: Hertzsprung-Russeh (HR) and G-M diagrams - Magellanic Glouds 
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1 INTRODUCTION 

Owing to its proximity, significant mass, and ongoing star 
formation, the Small Magellanic Cloud (SMC) represents a 
fundamental laboratory for studies of stellar and interstel¬ 
lar medium (ISM) processes. It contains a rich population 
of star clusters, associations, stellar pulsators, primary dis¬ 
tance indicators, and stars in short-lived evolutionary stages 
( Westerlund|1990 Bica Schmitt 1 1995 Bolatto et al.|2007 


SMC in terms of SFH and geometry. It constitutes the nat¬ 


ural follow-up of the papers by Kerber et al. (2009) and 


Soszyhski et al.|2010a|b 2011| ), with metallicities systemat¬ 
ically lower than those observed in the Large Magellanic 
Cloud (LMC) and Milky Way (MW) galaxies. 

The SMC is probably the largest among the Local 
Group galaxies in which a disk cannot be identified at first 
sight. Its optical appearance on the sky is dominated by an 
elongated bar-like structure in the direction NE-SW, from 
which a prominent Wing departs to the east, in a direc¬ 


tion that joins the Magellanic Bridge and the LMC (Nidever 
et al.|[20TT ). In addition, there have been consistent claims 
of depth structures being detected along different lines-of- 


sight ( Gardiner HawkinslpTM^ Nidever et al.|[20T3 ), and 
the growing awareness that the Bridge is, at least partially, 
a tidally-stripped component extending towards the LMC 
( Harris|20Q7 Nidever et al.|2013 ). Needless to say, any work 
on stellar or ISM physics can only be complicated by this 
particular morphology. 

The star formation history (SFH) of the SMC also car¬ 
ries along the signs of a disturbed past, probably marked by 
the interactions with the LMCQ The presence of periods of 
enhanced star and cluster formation has for long been rec¬ 
ognized, as well as a relative paucity of very old populations 


(e.g. Harris & Zaritsky|2004| 

de Grijs & Goodwin|2008 Noel 

et al.|2009| 

Piatti||2011c||Cignoni et al.||2012[ |2013a||Weisz 

et al.j 2013 

^ This disturbed past became even more evi- 


dent with the discovery that a fraction of SMC stars has 
been captured by the LMC ( Olsen et al.|20lf ). 

The VISTA Survey of the Magellanic Clouds (VMC; 
Cioni et al. 2011) was designed to greatly enhance our knowl¬ 


edge of the Magellanic system, via near-infrared photometry 
of the bulk of their stars, down to depths reaching the oldest 
turn-offs in both the LMC and SMC ( Kerber et al.|[2009 ). 
The advantage over previous deep-and-wide surveys resides 
mainly in the use of near-infrared wavelengths, that provide 
a view of the stellar populations almost unaffected by dust. 
The near-infrared light is also dominated by cool giants and 
subgiants, hence better sampling the intermediate-age pop¬ 
ulations over the young ones. In addition to depth, VMC 
provides twelve-epochs photometry in the Kg-band, allowing 
us to derive tighter period-luminosity relations for variables 


and more stringent constraints on the distances (e.g. Ripepi 
et al. 2012) and possible structures along the line of sight 
( Moretti et al.|2014 ). 

In this paper, we analyse the VMC photometry of the 


Rubele et al. (2012), in which the method for SFH-recovery 
from VMC data was first planned and tested, and the first 
analysis of the LMC data was performed, respectively. Since 
then, an impressive amount of data has been collected for 
the SMC galaxy, allowing us to perform a similar analysis. 
Present data for the SMC comprise most of its stellar mass, 
all of its bar, and a significant fraction of its Wing, as well 
as a few outer regions. 

The structure of this paper is as follows. Sect.|2jpresents 
the new VMC data and the PSF photometry. Sect7[^ details 
the method adopted for the derivation of the SFH and AMR, 
and their errors, as well as the best-fitting mean distances 
and reddenings for every subregion of a tile. Sect. [^discusses 
all the derived trends, also in comparison with previous re¬ 
sults in the literature. Finally Sect. [^summarizes the main 
results and draws a few concluding comments. 


2 DATA AND PHOTOMETRY 
2.1 Selected VMC data 

All data used in this work come from the VMC survey. 


which is extensively described in Cioni et al. (2011). We 


use the vl.l and vl.2 VMC data retrieved from the VISTA 
Science Archive (VSA; Hambly et al. ] |2004| p| More specifi¬ 
cally, we start from the pawprints already processed by the 
VISTA Data Flow System (VDFS; [Emerson et al.||2004 ) in 
its pipeline (Irwin et al. 2004[ ). Individual pawprints are com¬ 
bined into deep tiles as described in Appendix [A] The pho¬ 
tometric zeropoints are discussed in Appendix [B] 

In this work we deal with the the ten SMC tiles for 
which we already have collected a large fraction/the whole 
set of Kq multi-epoch photometry. They are illustrated in 
Fig.[^and listed in Table[^ The background image in Fig.[^ 
is a density map of all VMC sources with Ag < 18 mag 
and Kq errors smaller than 0.2 mag. Since this magnitude 
cut includes the red clump (RC) and the upper part of 
the red giant branch (RGB), the map is dominated by the 
intermediate-age and old stellar populations. It is immedi¬ 
ately evident that our ten tiles cover most of the SMC main 
body, including all of its dense bar, which appears in Fig.[^ 
as a bean-like blue region across the tiles SMC 3_3, 4_3, 4_4, 
5-3 and 5_4. In addition, the region analysed includes a large 
fraction of the so-called Wing (tiles SMC 3_5, 4_5 and 5_6), 
and more external regions towards the N and NE of the SMC 
(tiles SMC 6-3 and 6_5). We can confidently state that we 
have observed the majority of the SMC stellar mass, and the 
bulk of the young star formation. As we will see later, these 
observations also span a wide variety of SMC environments 
and distances. 


^ Although possible past interactions with the Milky Way are 
presently disfavoured ( |Besla et al.|[2007| , they might have oc¬ 
curred at ages larger than 4 Gyr, depending on the precise mass 
of the galaxies involved (|Kallivayalil et al.||2013| |G6mez et al.| 
2014]). 

^ i’Jie presence of a very old population is demonstrated by the 
presence of 2475 RR Lyrae ( [Soszyhski et al.]|2010b| ) detected by 
the OGLE III survey. 


2.2 PSF photometry and artificial star tests 

We have produced homogenised deep tile images and per¬ 
formed PSF photometry on them, as described in Ap- 
pendix[Aj Afterwards we correlated the three bands {Y JKs) 

^ http://horus.roe.ac.uk/vsa/ 
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Table 1. VMC tiles used in this work. 


Tile 

a (deg, J2000) 

5 (deg, J2000) 

Completion in RTg-band^ 

Comments 

SMC 3_3 

11.1995 

-74.2005 

92% 

S extreme of bar 

SMC 3_5 

21.9762 

-74.0018 

100% 

S part of Wing 

SMC 4_3 

11.2810 

-73.1094 

88% 

SW part of densest bar 

SMC 4_4 

16.3303 

-73.0876 

92% 

central bar, slightly towards Wing 

SMC 4_5 

21.2959 

-72.9339 

92% 

central part of Wing 

SMC 5_3 

11.2043 

-72.0267 

92% 

NW of densest bar 

SMC 5_4 

16.1918 

-71.9850 

100% 

NE part of densest bar 

SMC 5_6 

25.4401 

-71.5879 

88% 

~ 4 deg E of main body 

SMC 6_3 

11.4311 

-70.9266 

78% 

~ 1.5 deg NW of main body 

SMC 6 5 

20.4138 

-70.7601 

78% 

~ 3 deg NE of main body 


^ Here 100% completion corresponds to 12 epochs in the Xg-band or at least 9000 sec of integration time. 


-69^-30 

7(]"00 

-70^30 

7L"00 

-7r'30 

-72"00 

-72''3n 

-73"00 

- 73"30 

-74"00 

- 74^30 



Ihan lh50 lh40 IhlO Ihla LhLD LhOD OhW 0h4fl Dh30 {>h20 


Figure 1. Sky position of all tiles analysed in this work, as listed in Table ^ Each tile is subdivided into 12 small subregions, as 
illustrated here for the tile 5_6. The plot includes only the tiles with a high fraction of the planned VMC observations already completed; 
the complete, final coverage of the SMC by the VMC survey can be appreciated in |Cioni et al.| ( [2011| ). 


photometry using a 1^^ matching radius generating a multi¬ 
band catalog. Finally, we corrected our photometry for the 
aperture using as reference the VS A data release vl.2 (see 
Cross et al. poTl] and [Irwin et al.|2004| for details). 


Figure shows the CMDs in Y — Ks vs. iCg for the ten 
SMC tiles investigated to date. The different colours map the 
density of stars in the CMD. These CMDs will be further 
discussed in the next sections, together with the derived 
SFHs. 


Very extensive artificial star tests (ASTs) were per¬ 
formed on the homogenized images, so as to completely map 
the distributions of photometric errors and completeness, as 
a function of colour, magnitude, and position. The process is 
the same as extensively described and illustrated in [Rubel^ 


et al. (2012), with the only difference being that the photom¬ 
etry was made on the entire tile (see Appendix 0- Suffice 
it to mention here that our typical 50% completeness limit 
in the outermost tiles turns out to be at about 22.3, 22.1 
and 20.8 mag in the V, J and iCg filters, respectively. In 
the more central tile SMC 4_3, the same limits are found at 
about 21.3, 21.1 and 20.6 mag, respectively. 


In the following, we will only use VMC data and CMDs 
at magnitudes brighter than Ks < 20.5 mag, hence ensuring 
that completeness never falls below the 50% level. 
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Figure 2. PSF photometry for the ten tiles analysed in this work, in the Ks versus Y — Ks diagram. Colours from dark blue to red 
code the increasing density of stars in the CMD while the black bars on the left-hand side show the zblcr errors in iFg-band. The grey 
dashed lines mark the 50 % completeness level in the Ks band. Finally, the dark dot-dashed lines show the main sequence to RGB 
section of a 10 Gyr old isochrone of metallicity [Fe/H] = —1, shifted by a distance modulus of (m —M)o = 18.92 mag and an extinction 
of Ay = 0.35 mag, for comparison. 


2.3 Defining subregions 

All tiles, covering total areas of ~ 1.5 deg^ each, were sub¬ 
divided into twelve subregions of 21.0'x 21.5' (~0.12 deg^) 
for the SFH analysis. They are numbered from G1 to G12 as 
already illustrated in Fig. Gentral coordinates are listed 
in Table [ 3 ] below. G1 is located at the SE extreme of each 
tile, G4 at the SW, G9 at the NE, and G12 at the NW. 

The pawprints contributing to the corner subregion G1 
(see Eig. include a contribution from the “top” half of 


the VIRGAM detector number 16, which is known to show 
a significantly worse signal-to-noise ratio than the other de¬ 
tectors (see Rubele et al. 2012 for details). The effect is 
negligible in iCg, small in H and becomes more noticeable in 
the bluest VISTA passbands, i.e. in J, Y and Z. Therefore 
part of subregions G1 of all tiles were excluded from our 
analysis. 
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Figure 3. From left to right: Hess diagrams for the data, best 
fitting model, and chi-square, for the subregion G8 of tile SMC 
5-4. The top panels are for the J — Ks vs.Ks CMD, the bottom 
ones for Y — Ks vs. Kg CMD. In the left and central panels, the 
colour scale indicates the number of stars per bin. 


3 DERIVING THE SFH 
3.1 Overview of the method 


The method for deriving the SFH from VMC data has been 


already described extensively in Rubele et al. (2012). As a 


short summary, we recall that, for every subregion of a tile: 

(i) We start assuming an absolute distance modulus, 
(m —M)o and V-band extinction. Ay- 

(ii) “Partial models” are derived for the entire range of 
metallicities and ages of relevance (see next section). Partial 
models are simply synthetic stellar populations each cov¬ 
ering a small bin of age and metallicity, shifted to the de¬ 
sired distance and extinction, and representing a given total 
mass of formed stars. The extinction in each VISTA pass- 
band is derived from the assumed Ay using the extinction 
coefficients computed from the |Cardeffi, Clayton fc Mathis 
(1989) extinction curve with Ry = 3.1 (see Girardi et al. 
2008 for details), namely Ay = 0.385Aa/, Aj = 0.2S3Ay 
and Aks = 0.114Av|^. 


^ We choose to express the extinction in the V band just for 


(iii) Finally, partial models are “degraded” to the con¬ 
ditions of the actual observations by applying the distri¬ 
butions of photometric errors and completeness implied by 
the ASTs, and translated into Hess diagrams, i.e. maps of 
the stellar density across the Ks versus J — Ks and Ks ver¬ 
sus Y — Ks CMDs. This involves a few steps: First, high- 
resolution Hess diagrams, with a resolution of 0.05 mag 
in magnitude and 0.04 mag in colour, are built from the 
synthetic stellar populations. For each one of these small 
boxes, all artificial stars of similar colour and magnitude 
are looked for in the AST database, in order to build two- 
dimensional histograms of the differences between the input 
and output magnitudes and colours, and compute the com¬ 
pleteness fraction per box. The high-resolution Hess diagram 
is then degraded by applying these error distributions and 
completeness on a box-by-box basis, which has the effect of 
both spreading stars into the neigbouring boxes in the CMD, 
and reducing their numbers. Finally, the high-resolution de¬ 
graded partial models are converted into Hess diagrams of 
normal resolution (namely 0.15 mag per 0.12 mag), which 
are used in the subsequent steps of the analysis. 

(iv) The linear combination of partial models that best-fit 
the observed Hess diagrams is found by the StarFISH opti¬ 
mization code by Harris V Zaritsky (2001), suitably adapted 
to our case. An example of a best-fit model is presented in 
Fig. in The coefficients of this linear combination of par¬ 
tial models are straightforwardly converted into the SFH. 


mesures the residuals between the best-fit model and the 
data. 

(v) The same process is repeated over a range of different 
(m — M)o and Ay values. This allows us to pin down the 
best-fit SFH, distance and extinction for each subregion, by 
simply looking for the minimum Xmin- 

(vi) One hundred synthetic realizations of the best-fit 
model are generated and are analysed via the same method. 
The dispersion in the resulting coefficients provides us the 
confidence levels of our best-fit SFH, {m — M)o and Ay. 


Several aspects of this process will be either better detailed 
or simply illustrated in the following, whenever relevant to 
understand the particular results we find for the SMC. 


3.2 Partial models 


Partial models for this work have been derived from PAR¬ 
SEC vl.O evolutionary tracks ( [Bressan et al. 2012), which 
represent a major revision and update of the previous sets 
of Padova evolutionary tracks (namely Marigo et al.||2Q08 


Girardi et al.||2010 ) used in [Rubele et ai.| j2012). The most 


relevant changes, in the context of this work, are 


the adoption of Caffau et al. (2011 and references 


convenience in the subsequent comparisons with other authors 
(Sect. |4^ , although we will be actually dealing with the extinc¬ 
tion in near-infrared pass-bands. Notice that the average LMC 
and SMC extinction curves are essentially indistinguishable from 


[don et al.||20Q3| figure 10), so that the choice of a near-infrared 
reference wavelength for the extinction, instead of Ay, would lead 
essentially to the same results. 


the MW one for all wavelengths redward of ~ 5000 A (see e.g. 


Cor- 
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Table 2. Grid of SMC stellar partial models used in the SFH 
recovery. 


log(t/yr) 

[Fe/H]i 

[Fe/H] 2 

[Fe/H],, 

[Fe/H]4 

[Fe/H]., 

6.9 

-0.10 

-0.25 

-0.40 

-0.55 

-0.70 

7.4 

-0.10 

-0.25 

-0.40 

-0.55 

-0.70 

7.8 

-0.10 

-0.25 

-0.40 

-0.55 

-0.70 

8.1 

-0.10 

-0.25 

-0.40 

-0.55 

-0.70 

8.3 

-0.20 

-0.35 

-0.50 

-0.65 

-0.80 

8.5 

-0.20 

-0.35 

-0.50 

-0.65 

-0.80 

8.7 

-0.20 

-0.35 

-0.50 

-0.65 

-0.80 

8.9 

-0.40 

-0.55 

-0.70 

-0.85 

-1.00 

9.1 

-0.55 

-0.70 

-0.85 

-1.00 

-1.15 

9.3 

-0.55 

-0.70 

-0.85 

-1.00 

-1.15 

9.5 

-0.70 

-0.85 

-1.00 

-1.15 

-1.30 

9.7 

-0.85 

-1.00 

-1.15 

-1.30 

-1.45 

9.9 

-1.15 

-1.30 

-1.45 

-1.60 

-1.75 

10.075 

-1.45 

-1.60 

-1.75 

-1.90 

-2.05 


therein) solar composition, with a present-day solar metal 
content oi Z — 0.0152; 

• the recalibration of the Solar Model so as to meet the 
stringent constraints imposed by helioseismology and by the 
radiometric age of the Solar System; this implies allowing 
microscopic diffusion to operate in low mass stars; 

• revised limits for the mass interval for which overshoot¬ 
ing increases its efficiency, in low-mass stars. 


The main consequence of these novelties are modest (though 
systematic) changes in the age and metallicity scale of stellar 
populations, at ages exceeding a few Gyr. The age scale is 
presently being tested on LMC/SMC intermediate-age star 
clusters for which we have HST data of excellent quality 
e.g. iGirardi et ah 2013[ [Gorrenti et al. 2Q14| Goudfrooij 


et al.||2014), so that we postpone further comments on this 


topic to future papers. Also, we recall that the present re¬ 
lease of PARSEG tracks (vl.l as of this writing) already 
includes modifications in the prescriptions for microscopic 
diffusion in low-mass stars, with respect to vl.O. We have 
checked that these changes are simply negligible in the range 
of magnitudes and colours sampled by the VMG data for the 
SMG. 

Table displays the location of our partial models on 
the [Fe/H] vs. log(t/yr) plane. Age bins are defined as in 
Rubele et al.| ( |20l'^ ; they are, in general, at equally spaced 


intervals of Alog(t/yr) = 0.2 dex, so as to reflect the fact 
that our age resolution is better at young ages. For the very 
youngest age bins, however, wider age intervals are used, so 
as to increase the statistics in places where the star counts 
are naturally very small. 

A total of 5 partial models are created for every age bin. 
They have metallicities covering an interval of A [Fe/H] = 
±0.3 dex with respect to the reference metallicity. The lat¬ 
ter is derived from the Piatti (2011c) AMR, which results 
from observations of SMG star clusters. Metallicity bins are 
separated by 0.15 dex from each other. It is important to 
remark that the Piatti (2011c) AMR is used just to define 
the location of partial models, but it has virtually no effect 
in determining the final AMR which will be derived from 
our data (see Sect. |4.6| later). 

To this set, we add one partial model representing the 
Milky Way foreground. This is generated by running the 


TRILEGAL Milky Way model ( |Girardi et al.|2005|p 


for the 

central coordinates of each subregion, and for the total area 
it represents. 

We also note that partial models are computed for the 
entire range of distance moduli and extinctions explored in 
this work, namely 18.6 < (m —M)o/mag < 19.3 and 0 < 
Av/niag<0.8. At the high end of this distance/extinction 
ranges, the 10 -Gyr main sequence turn-off might be located 
about 0.5 mag fainter than the Ks < 20.5 mag limit being 
adopted in the derivation of the SFH. Anyway, this magni¬ 
tude interval includes a good fraction of the subgiant branch 
for all ages, hence ensuring a good sensitivity to the old SFH. 
Moreover, we recall that the constraints on the oldest SMG 
populations do not depend on the main sequence and sub¬ 
giant branches only, but also on their contribution to the 
RGB, RG, and horizontal branch, as demonstrated by |Dol-| 


phin ( 2002 ), and tested in the preparatory work for VMG 


by |Kerber et al.| ( |2009 ). 

Finally, it is worth mentioning that all models used in 
this work were produced using the Vegamag definition of 
photometric zeropoints, according to the procedure detailed 
in Girardi et al. ( 2002 ). We seek for a good consistency be¬ 


tween this definition and the zeropoints derived during the 
calibration of VMG data, as discussed in Appendix For 
the moment, we just warn the reader that there might exist 
residual offsets between these two realisations of VISTA ze¬ 
ropoints, that might primarily appear as small offsets in the 
values of distances and extinctions derived in this work. 


3.3 Overview of the results for the SFR(t) and 

[Fe/H](i) 

The next steps leading to the final best-fitting SFHs are 
illustrated in Fig. which refers to the analysis of subre¬ 
gion G 8 of tile SMG 5_4. Essentially, we run the StarFISH 
SFH-recovery software over a wide-enough grid of distance 
modulus and extinction values, (m —M)o and Ay, as illus¬ 
trated in the left panel. For each point in this grid, the code 
returns not only the best-fitting coefficients for the SFH, but 
also the Xmin- Us minimum value allows us to identify the 
best-fitting {{m — M)o, Ay) pair, which in this case falls at 
{m — M)o = 18.89 mag and Ay = 0.53 mag. The coefficients 
of this best-fitting solution contain all the information nec¬ 
essary to express the SFH in the age-metallicity plane. For 
a more convenient representation of the results, for each age 
bin we derive the summed coefficients, producing the star 
formation rate function SFR(t) in units of M 0 yr“^, and 
the average metallicity, hence deriving the age-met allicity 
relation, [Fe/H](t). They are illustrated in the right panel of 
Fig. ID 

We also evaluate the random and systematic errors 
in the SFR(t) and [Fe/H](t) functions. One hundred sim¬ 
ulations are produced using the TRILEGAL tool and the 
best-fitting solution, and then degraded using the ASTs 
and re-analysed by the SFH-recovery software. The scat¬ 
ter among these 100 solutions gives the SFH random errors. 
The distribution of the Xmin among these simulations al¬ 
lows us to identify the Xmin values that correspond to the 
la and 3a (68 and 99.7 %, respectively) confidence levels. 


^ http://stev.oapd.inaf.it/trilegal 
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Figure 4. The SFH results for subregion G8 of tile SMC 5_4, which corresponds to a relatively dense region of the SMC bar. The left 
panel shows the map of Xmin derived from StarFISH over a wide range of distance moduli and extinctions. The continous line delimits 
the 3(7 (99.7%) confidence level area of the best-fitting solution, while the dashed line delimits the Icr (68%) confidence level. The right 
panels shows the best fitting solution in the form of the SFR in units of Mq yr“^ as a function of the logarithm of age (blue histogram, 
top panel), for all age bins represented in the partial models, together with the averaged metallicity (black dots, middle panel). In the 
middle panel, the small green dots mark the central [Fe/H] and log(t/yr) values for which partial models are defined. In both panels, 
the vertical error bars mark the random errors in the SFR(t) and [Fe/H](t) relations, while the grey-shaded areas delimit the systematic 
errors. 


Table 3. Mean distances and extinctions for each subregion. A complete version of this table is provided in electronic form on the journal 
website. 


Tile 

SMC 

subre¬ 
gion C 

a 

(deg) 

s 

(deg) 

Mo 

(mag) 

^IJ,Q^loW 

(mag) 

(mag) 

d 

(kpc) 

^d,low 

(kpc) 

^d,up 

(kpc) 

(mag) 

^ Ay ,1ow 

(mag) 

^Ay ,up 

(mag) 

3_3 

1 

13.18529695 

-74.59449494 

18.92 

0.10 

0.08 

60.81 

2.68 

2.81 

0.426 

0.051 

0.099 

3_3 

2 

11.94330634 

-74.61137246 

18.90 

0.07 

0.08 

60.26 

2.01 

2.08 

0.454 

0.079 

0.046 

3_3 

3 

10.59062191 

-74.60080543 

18.97 

0.10 

0.08 

62.23 

2.78 

2.91 

0.382 

0.057 

0.043 

3_3 

4 

9.23135534 

-74.58694737 

19.02 

0.07 

0.08 

63.68 

2.09 

2.16 

0.369 

0.044 

0.031 

3_3 

5 

13.27320506 

-74.20450581 

18.93 

0.08 

0.07 

61.09 

2.11 

2.18 

0.409 

0.035 

0.040 

3 3 

6 

11.95494707 

-74.20545383 

18.98 

0.05 

0.07 

62.52 

1.54 

1.58 

0.416 

0.041 

0.084 


Systematic errors instead are evaluated by simply examin¬ 
ing the total range of all the SFH solutions inside the la 
area of the ((m — M)o, Ay) plane. The detailed procedure 
has been more extensively described in Kerber et al. (2009) 
and Rubele et al.| ( |2012| ) . 


Results for the 120 subregions analyzed in this paper are 
presented in Figs. and In these figures, we present the 
SFR(t) and [Fe/H](t) resulting for all tiles and subregions, 
displaying them in tabular form, in the same way as they 
appear in the sky. The figures illustrate clearly the strong 
correlation between the SFH results for nearby subregions 
and tiles, especially for ages larger than log(t/yr) > 9. The 
continuity in the SFH between neighbouring tiles is particu¬ 
larly encouraging, since - contrary to subregions of the same 
tile - they do not necessarily share the same photometric 
quality and completion in the iFg-band (see Table [^. 

Another feature evident from Fig. I^is that the SFR(t) 
for the outermost SMC regions is overall quite similar, as 
can be appreciated comparing panels in the extreme corners 
of the surveyed area. This “external” SFR(t) is composed of 


a main period of star formation between log(t/yr) = 9.2 and 
10.0, plus a trace of star formation at young ages. In most 
of the external subregions, however, the young SFR is not 
significant compared to the error bars, so that we have only 
upper limits to it. Fields across the main SMC bar and on 
the Wing, instead, have a much more complex SFR(t), with 
seemingly different bursts of star formation appearing along 
their histories, including very significant events of young star 
formation. Moreover, we can see that what appears as a sin¬ 
gle period of steady SFR for all ages log(t/yr) > 9.2 in 
the extreme NW and SW of the survey area, in the East¬ 
ern most remote areas appears more like a double-peaked 
SFR(t), with maxima at log(t/yr) = 9.5 and 9.9. These two 
peaks are more evident in the NE tile SMC 6_5. These fea¬ 
tures will be discussed further in Sect. 14.51 below. 
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Figure 5. An overview of the SFR(t) for all subregions analysed in this work. The small panels present both the best-fit SFR(t) (red 
line), and the Icr confidence level (shaded grey area) for every subregion. They are displayed in tabular form, in the same way as they 
appear in the sky (see Fig. Q. The axis limits are not presented in every small panel, because this would make the plot too crowded; 
anyway the axis limits can be easily derived with the help of the enlarged panel at the bottom-left, which serves as an example: All 
panels have the same limits in the ordinate, i.e. 6.2 < log(t/yr) < 10.5. The limits for the SFR in the abscissa wildly vary between 
panels: they are set between 0 and a number which is presented inside each panel, with two significant digits; this number is in units of 
lO-'^Moyr-i. 


4 DISCUSSION 

Although we obtain best-fit solutions and confidence inter¬ 
vals simultaneously for all parameters - SFR(t), [Fe/H](t), 
(m —M)o and Ay - it is convenient to discuss them sepa¬ 
rately in the following, postponing the general discussion of 
all trends to the final subsection. 


4.1 The SMC geometry 


As already anticipated in the introduction (see also Sect. 4.2 
below), several authours have derived significant depths 
along the SMC lines-of-sight. We, instead, have used the 
standard approach of assuming a null depth in the SFH 
analysis. This approximation might have two different im¬ 
plications. First, one may wonder whether the best-fitting 
distance and SFH we derive for each subregion are represen¬ 
tative of the mean distance and SFH that would have been 
derived if a significant depth was adopted. This has been 


already verified by Harris & Zaritsky (2004), who simulated 


populations with depths up to 12 kpc (0.2 mag) and then 


analysed them with StarFISH without accounting for the 
spread in distance. They conclude that the SFH solutions 
were the same within the errors, and moreover that there 
were no indications of a significant depth in their results for 
the SMC. Since we are using the same software and method, 
we essentially endorse their conclusions and start assuming 
no depth a priori. The second implication of our null-depth 
approximation is that we will be constrained to analyse the 
SMC geometry more in terms of a “thin surface projected on 
the sky” rather than in terms of a complex 3D distribution. 
Changing this approach requires a dramatic change in the 
methods used in the SFH analysis, which will be pursued in 
subsequent works. 

That said, the best-fitting and 68% confidence intervals 
of the distances and extinctions are presented in Table 
(a complete table is provided in electronic format on the 
journal website). 

All distances are plotted as a function of equatorial co¬ 
ordinates in Figs. (left panel) and which indicate sig¬ 
nificant variations in its mean distance across the SMC and 
especially along the E-W direction. While central parts of 
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Figure 6. The same as Fig. but for the AMR, [Fe/H](t). This time, all panels present the same limits, as illustrated in the enlarged 
panel at the bottom-left. 
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Figure 7. Maps showing the best-fit distances (left panel) and the widths of its 68% confidence intervals (right panel) for each subregion. 


the galaxy are found at distances between 58 and 60 kpc, 
the tiles SMC 3_3 and 6_3 are clearly farther out at ~62 kpc, 
while the easternmost tile SMC 5_6 - the one closest in the 
sky to the LMC, and possibly associated with the Bridge - 
is found at a distance of ~54 kpc, which is more comparable 
to that of the LMC rather than the SMC. As can be appre¬ 


ciated in Fig.|^ there is a clear trend of increasing distances 
in the E-W direction. 

We translate this E-W trend of the distances into a 
simple disk model as in the case of the LMC (see [Rubele 


et al. 2012). This was performed for two different choices 


of the SMC centre: the apparent kinematic centre from 
Stanimirovic, Staveley-Smith & Jones ( 2004| S04), with 
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10 Rubele et al. 



Figure 8. Distances and their 68% confidence intervals as a function of coordinates a and 8. It is evident that the distances increase 
systematically with a, whereas there is no clear trend with 8. 


(ac = 16.25°, (5c = —72.42°), and the stellar density cen¬ 
tre indicated by the VMC data. The latter corresponds to 
(ac = 12.60°, (5c = —73.09°), and was determined as illus¬ 
trated in Fig. We remark that it almost coincides with 
the centre of the K and M giants determined by |Gonidakis| 
et al. ( 2009| ), with (ac = 12.75°,(5c = —73.10°). A weighted 


least-squares fit of a plane to the distance data then provides 
a position angle of the line of nodes of Oq = 179.5 ± 1.3° and 
an inclination angle of z = 39.7 ± 5.5° in the case of the 
S04 centr(|3 and Oo = 179.3 ± 2.1° and i = 39.3 ± 5.5° in 
the case of our centre. The distances to the kinematic and 
stellar-density centres’ intersection with the best-fit plane 
are (m — M)o = 18.87 ± 0.01 mag (59.48 ± 0.28 kpc) and 
(m — M)o = 18.91±0.02 mag (60.45±0.47 kpc), respectively. 
These two cases are illustrated in Fig.[^ Residual distances 
from these best-fit planes are reasonably high, amounting to 
up to ~ 2 kpc in both directions. However, it is only in the 
easternmost regions that the approximation of a disk clearly 
fails, with residuals approaching ~ —3 kpc. Moreover, it can 
be noticed that at about 1.5 kpc to the east of the SMC stel¬ 
lar centre (which coincides with the kinematic centre) the 
best-fit plane is located systematically in front of the data 
points by about 1 kpc. These deviations are generally larger 
than the la errors in the distance determinations, as can 
be appreciated in the bottom panels of Fig.[^ The general 
picture we get is that the SMC is very distorted, and if its 
structure is to be approximated by a disk, it is a significantly 
warped one. Another interesting aspect is that the position 
angle of the line of nodes of the best-fit plane appears not 
to be related with the major axis along which most of the 
SMC young stars and gas are distributed (i.e. the so-called 
SMC bar), which has a position angle of ~ 40°, and which 


also coincides with the kinematic position angle (cf. Stan- 
imirovic, Staveley-Smith Jones|2Q04 ). Similar results have 


been obtained recently by 


Tatton 


(2014), from the analysis 


of the RC magnitude in VMC data. 

It should be stressed that the choice of fitting a disk is 


® ^0 and i are defined as in |van der Marel| ( |20Ql| , i.e. is mea¬ 
sured starting from the north in counter-clockwise direction, and 
the near side of the SMC is at — 90°. 



_._i_■_I_■_I_■_k_■_I_._I_■_ i _._I_■_ i _■_I_■_ I 

L7j 6 16.14 L6 L52 14.4 L3.6 12.14 12 LL.2 10.4 y.6 14.14 

RA 


Figure 9. Stellar density map of the central SMC, where stars 
are selected as described in Sect.^ The black contour line is the 
density limit inside which the density variation is less than 10% 
away from the peak value; this region is used to estimate the 
SMC stellar-density centre (pink cross), located at (a = 12.60°, 
8 = —73.09° and (m — M)o = 18.856 mag (d = 59.04 kpc). 


not determined by a belief that the SMC has (or is) a disk; 
it is mainly driven by two facts: (1) it is the simplest surface 
that can be fit to a table of distances as a function of coor¬ 
dinates, as the one which we have derived; and (2) because 
finding the mean disk parameters allows us to easily compare 
our distance results with a series of other similar derivations 
of the SMC “best-fit plane” available in the literature. The 
two main parameters we find, namely the inclination and 
position angle, could also be interpreted in terms of other 
simple 3D configurations like, for instance, oblate or prolate 
spheroids. 

Regarding the comparison with other results in the liter¬ 
ature, it is generally found that both the old and intermedi¬ 
ate age populations (ages > 2 Gyr) are distributed in a pres¬ 
sure supported, spheroidal/ellipsoidal distribution ( jSubra- 


manian Subramaniam|2012 and references therein). The 


young stars, like classical Cepheids, are found to be in a 
disk. In the present work we do not use any specific distance 
tracer, but the use of near-infrared data implies that our 
distances are generally more weighted by the intermediate- 
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distance (kpc) along the line of maximum gradient (~ E-W) 



distance (kpc) along the line of maximum gradient (~ E-W) 


Figure 10. The two upper panels show the distance of the subregions as projected perpendicularly to the line of nodes, i.e. along the line 
of maximum distance gradient, which is close to the E-W direction. This line (dashed) is determined from fitting a plane to the distance 
data, for two different choices of the SMC centre (marked by a blue diamond): adopting either the kinematic (S04, left panels) or stellar 
density centre (this work, right panel). The dotted lines represent a projected disk thickness of 1 kpc (or 0.82 kpc linear distance) to 
both sides of the best-fitting plane. The bottom panels show the residuals with respect to best-fit plane, expressed in units of the Icr 
uncertainty in the distance determination. 


age populations than by young stars. So, it is initially a 
surprise that we find indications of a distorted disk, rather 
than the more smooth features that would be expected from 
a spheroidal/ellipsoidal distribution. Moreover, it also sur¬ 
prising that the structure we (and also Tatton 2014) ob¬ 
tain is very different from the one derived from the clas¬ 
sical Cepheids, which indicate a much more inclined disk 
with z ~ 65 — 70° (Caldwell & Coulson 1986). Cepheids 
also possess a very different position angle of line of nodes, 
PA ~ 145 — 150° (see also Laney Stobiep986 Haschke, 
Grebel Duffau]|2012| ). 


4.2 Possible depth structures 

The total distance intervals corresponding to the 68% con¬ 
fidence level of the best-fit distances are illustrated in the 
right panel in Fig. One can see clearly that such widths 
are larger towards the SE side of the SMC, and especially 
at the SE border of tile SMC 3_5, which has a distance con¬ 
fidence interval about 8 kpc wide. Most of the other fields 
have best-fit distances constrained to within 5 kpc - with a 
few fields having distances constrained to be inside < 1 kpc. 
One can also notice that there is a mild trend of the total 
distance interval increasing gradually towards the SE of the 
SMC. 

The interpretation of this increased distance confidence 
level intervals is not straightforward, but two possible causes 
for them are: (1) indications that these areas present large 
distance spreads along the line of sight, or (2) that the 
smaller stellar density towards these outer regions are caus¬ 
ing sitematically poorer fits of the SFHs and distances. In¬ 


deed, IKerber et al.|[200^ have shown that the SFH errors 
are primarily a function of stellar density and completeness. 
However, it is striking that the tiles SMC 5_6 and 6_5, to¬ 
wards the NE, present similar stellar densities (see Fig, 
and smaller completion than tile SMC 3_5 (see Table [y , 
but no indication of an increased distance confidence level 
interval. This lead us to suspect that the SE areas of the 
SMC might indeed present significant depths along the line 
of sight. This suspicion cannot be confirmed at the moment, 
since our method does not assume a distance dispersion. A 
subsequent work will be dedicated to this point. 


In this context, however, it is worth recalling that many 
other authours have found indicated distance variations 
and/or a sustantial depth across different lines-of-sight of the 
SMC, starting with the early work by [Gardiner Hawkins] 
( 1991| ). The literature of the topic is highly varied and based 
on very heterogeneous data. The latest works on the subject, 
however, tend to use the RC position and widths as provided 
by major surveys such as the Magellanic Cloud Photometric 
Survey (MCPS) and the Optical Gravitational Lensing Ex¬ 
periment (OGLE), as the primary probe of the SMC geom¬ 
etry. Depths derived in this way often reach values between 
10 and 23 kpc ([Gardiner Hawkins||1991 Subramanian & 


Subramaniam[|2009[ Nidever et al.||2013 Tatton|2Q14 ). 

Interestingly, using the position and width of the RC. 


Nidever et al. (2013) find that the eastern side of the SMC 


has a bimodal distribution of distances, with a component 
located at just 55 kpc, far closer than the main SMC body 
which they locate at 67 kpc. Our best fit-distances (left panel 
in Fig.|^ essentially confirm the presence of this close com¬ 
ponent of the SMC, which is clearly seen in tile SMC 5_6 
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12 Rubele et al. 


(Fig. 13 left panel) for typical distances of 56 kpc. How¬ 
ever, our distance maps do not produce evidence for the 
farther component (expected at > 60 kpc) at this position. 
This could simply be caused by the nearby component being 
denser than the farther one, so being the only one being sin¬ 
gled out by our method. Indeed, examination of the CMDs 
for the SMC tile 5_6 (Fig. reveals the presence of an ex¬ 
tended RC, with the faintest half being less populated than 
the brightest; this latter probably corresponds to |Nidever| 


et al. (2013)’s distant component. 


Finally, before closing the discussion related to the 
distances, it is worth recalling that the distances deter¬ 
mined with our method are average values obtained assum¬ 
ing that in any subregion stars of all ages are at the same 
mean distances. The most representative ages sampled by 
this method vary a lot across the SMC, as we will see in 
Sect. |4.4] but it is fair to state that they are most sensitive to 
the intermediate-age populations that (if present) generate 
prominent RCs and RGBs in our near-infrared CMDs. How¬ 
ever, there are clear indications from other authors that pop¬ 
ulations of different ages have different geometries. For in¬ 
stance, based on the large sample of variables from OGLE HI 


Haschke, Grebel & Duffau (2012) find that the young clas¬ 


sical Cepheids are distributed with a position angle of the 
major axis ^ = 66 ± 15° and inclination z = 74 ± 9°, while 
the very old RR Lyrae variables are better described by a 
0 — 83 ± 21° and z = 7 ± 15°. Both distributions are clearly 
different from the mean geometry we hnd. This is not sur¬ 
prising for the young populations, which are concentrated 
on the SMC bar. The different results of the very old pop¬ 
ulation, instead, seem to indicate that this old population 
constitutes an extended halo (see also Nidever et al.|[20TT ) 
with little resemblance to the warped disk galaxy we derive. 

Needless to say, future analysis of the SFH based on 
VMC data will include a more detailed description of the 
distance distributions along the different lines-of-sights, and 
as a function of stellar ages. In this context, the present solu¬ 
tions constitute a valuable first-guess for this more complex 
analysis. Moreover, we are also studying the structure of the 
SMC as traced by classical Cepheids and RR Lyrae using the 
whole sample from OGLE HI plus EROS-2 variables which 
have a counterpart in VMC, and new period-luminosity re¬ 
lations derived from the multi-epoch VMC data (see e.g. 
Ripepi et al.||2012 2014 Moretti et al.||2014 ). This will 


able us to make a much more robust comparison between 
the geometry derived from the variables and from the SEH- 
recovery work. 


4.3 Large-scale extinction maps 

The extinction map derived from our best-fitting values of 
Av is presented in Eig. As can be appreciated, most of 
the external SMC regions are found at extinctions typically 
smaller than Ay = 0.5 mag, with minimum values being of 
Av — 0.35 mag. Only across the main bar, and in parts of 
the Wing regions in tile SMC 3_5, are the Ay values larger 
than this, reaching Ay > 0.65 mag for three subregions. 
The average extinction for the entire analysed area is of 
Ay = 0.47 ± 0.09 mag. We remark that the maximum value 
of Ay = 0.65 mag corresponds to extinctions of just 0.25, 
0.18 and 0.07 mag in the Liters V, J and iFg, respectively. 
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0 . 6 -^ 
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Figure 11. Map of best-fit extinctions Ay. 


and to colour excesses of 0.18 mag in Y — K and 0.11 mag 


in J — K (cf. item I (ii) I in Sect. 3.1). 

The comparison with extinction values from the lit¬ 
erature is not easy, given the very different methods and 
probes applied by different groups. The widely-used extinc¬ 
tion maps from Schlegel, Einkbeiner V Davis (1998), based 
on estimating dust column densities from dust emission, are 
not reliable for the central parts of the SMC; for its im¬ 
mediate surroundings, however, they indicate mean values 
close to Av^~0.12 mag, which probably represent the typi¬ 
cal values for the foreground MW extinction across the en¬ 
tire SMC. Similar dust maps were derived for the central 
SMC regions by Israel et al. (2010), using COBE/DIRBE 
and WMAP data; they provide a mean extinction of about 
Ay = 0.45 mag internal to the SMC. That means that their 
total average extinction values reach Ay ~ 0.57 mag, which 
compares very well with our results for the SMC subregions 
with the highest mean extinctions. Overall, however, |Israel| 


et al. (2010) extinctions appear larger than our values. This 


could be explained by the fact that these extinction maps 
are derived from emission-based estimates of the total dust 
column density, hence comprising the entire line-of-sight to 
the SMC galaxy, while our extinction values are a star-by¬ 
star weighted mean that reflect the mean values somewhere 
in the middle of the SMC stars. 

On the other hand, our average Ay values agree quite 
well with those derived by Zaritsky et al. (2002) from MCPS 
data, which span the range of 0.15 — 0.65 mag, when consid¬ 
ering both cool and hot stars. In their case, the extinction 
is probed using similar objects (the stars), but using a dif¬ 
ferent method, namely a star-by-star fitting of the spectral 
energy distribution. 


Haschke, Grebel & Duffau (2011) instead find an av¬ 


erage extinction of just Ay = 0.1 ± 0.15 mag (or Ey-i = 
0.04 ± 0.06 mag) adopting a theoretical mean unreddened 
colour for the RC stars, and Ay = 0.18 ± 0.15 mag (or 
Ey-i = 0.07 ± 0.06 mag) from RR Lyrae stars from the 
OGLE HI survey. A similar distribution of the extinction in 
the central part of the SMC was found in [Subramanian fc] 


Subramaniam (2012), also using the mean colour of the red 


clump. In both cases, the derived extinction values appear 
to be significantly smaller than ours. We should however 
keep in mind that these are mainly measuremements of the 
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internal SMC extinction, that is, the RC method is mostly 
sensitive to the variations of extinction with respect to a 
fixed reference value. Moreover, there is a significant differ¬ 
ence in assumptions with respect to the present work: in 


practice, both Haschke, Grebel & Duffau (2011) and Sub- 
ramanian & Subramaniam (2012 | ) assume that the RC has 


a constant intrinsic colour across the SMC. We, instead, 
assume that all CMD features (including the RC) vary as 
indicated by stellar models and following the age and metal- 
licity distributions of the subregions being sampled. Tracing 
back the origin of differences in the final Ay values, in these 
circumstances, is not straightforward. 

However, if we consider the two least-extincted north¬ 
eastern tiles (namely tiles SMC 5_6 and 6_5), there indeed 
appears to be a modest zero-point offset of ~ 0.2 mag in 
our Av values, compared to the canonical 0.12 mag derived 
from Schlegel, Finkbeiner & Davis (1998) for the SMC out¬ 
skirts. We note that such an offset, if real, could be eas¬ 
ily explained by offsets of the order of ~ 0.05 mag in the 
VMC photometric zeropoints, and/or in the stellar models 
employed here. Moreover, Peek V Graves (2010) show that 
differences in dust temperature may lead to variations in 
the Schlegel, Finkbeiner & Davis ( |1998 ) Eb-v values of a 
few hundredths of a magnitude (a tenth of a magnitude in 
Av)^ so that it is not even clear whether the problem is 
in our results or on the dust emission maps. Although this 
surely deserves further investigation, we recall that such er¬ 
rors would not affect other aspects of this work, since they 
would, prevalently, just cause a small systematic offset in 
our distance moduli (that is, a A(m —M)o ~ 0.05 mag). In¬ 
deed, all comparisons in the distances and reddenings we are 
making here are chiefly meant to explore how these quanti¬ 
ties vary across the SMC in a differential way. Future works 
by the VMC team will better address the problem of deter¬ 
mining the absolute values of distances and extinctions from 
these data. 



4.4 The look back time star formation rate and 
mass assembly tomography 

Figuresj^ and [^encode a lot of information about the SFH of 
the SMC, including their spatial distribution and derivations 
of the uncertainties. However, they do not provide a good 
visualisation of the relative spatial distribution of the dif¬ 
ferent populations; moreover, those are low-resolution maps 
whose pixel elements are equal to the subregions. In order to 
improve upon this point, we have transformed our SFH re¬ 
sults, together with the original photometry data, into high- 
resolution maps of the SFH and history of the SMC stellar 
mass assembly. 

Our stellar population maps are built as follows: 

(i) For each one of the 120 subregions analysed, we take 
the coefficients of the best-fitting solution and build the total 
Hess diagram of each age bin i and metallicity j. 

(ii) For each star in the photometric catalogue, we iden¬ 
tify its cell in the Hess diagram. The total number of stars 
in that cell is nobs- 

(iii) We assign to the star the relative probabilities of be¬ 
longing to each one of the 14 x 5 age-metallicity bins, by 


Figure 12. An intermediate step in the Stellar Populations CMD 
Reconstruction. The top-left panel shows the Hess diagram of 
the best-fit partial model of a subregion in tile SMC 3_3 with 
log(t/yr) = 10.075 and [Fe/H] = —1.45 dex. The top-right panel 
instead shows a simulated stellar population (dots) of exactly 
the same age and metallicity, but with the individual stars being 
coloured according to the probability of them belonging to that 
specific stellar population - as opposed to the probability of be¬ 
longing to all the other partial models that have been fit to this 
subregion. The bottom panels shows the same but for a young 
stellar population, with log(t/yr) = 8.7 and [Fe/H] = —0.50 dex, 
in the same subregion. 

simply computing the ratio 

^obs 

Pij — '^ij,bfm X 

2-^ij 

where nij,bfm is the number found in the best-fit model. 

(iv) We sum the pij of the five metallicity bins holding 
the same age z, hence deriving pi. 

(v) We randomly assign the star to an age bin, following 
the probability set by the distribution of pi . 

(vi) Finally, we build density maps containing all stars 
which have been observed. There is a map for each age bin. 

An intermediate step of this process is illustrated in 
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Figure 13. Detailed SFR(t) maps. Each panel shows the SFR intensity for a given age bin, derived as detailed in the main text. The 
mean log(t/yr) of each age bin is indicated above the top axis; they coincide with the values listed in Table 
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Figure |13| continued. 


Fig. |12[ which shows how a population of a given age and 
metallicity (that is, stars belonging to a single partial model) 
would appear in terms of pij. The top half of the plot shows 
how stars belonging to an old age appear, when translated 
into a probability of being truly old - or, more specifically, of 
belonging to the partial model with (i, j) = (14,1): they ap¬ 
pear very prominently (pi 4 ,i > 0.3) if they are at the faintest 
part of the red clump and horizontal branch, and moderately 
prominent (0.05 < pi 4 ,i < 0.2) over the RGB and reddest 
part of the turn-off, while stars at the bluest part of the turn¬ 
off turn out to be given a much lower probability of being 
truly that old. For a young population, instead, largest prob¬ 
abilities of being truly young are in general assigned to the 
stars at the brightest part of the He-burning sequence, and 
at the bluest part of the turn-off, as illustrated in the bot¬ 
tom panels of Fig. 12 for a population with (z, j) = (7,3). 


In our method, the same probability assignment is made to 
all stars observed by VMC. It is evident that the resulting 
population maps will heavily weight the stars whose ages, 
due to their particular position in the CMDs, have a high 
probability of not belonging to other age bins. 

The resulting spatial maps are presented in Fig.J^ for 
all 14 age bins we have defined. As the previous Fig.[5pthey 
also reflect the intensity of the best-fitting SFR(t), but with 
a significant advantage: across each subregion, this SFR(t) 
is spatially distributed in proportion to the number of stars 
that likely have that age, and the resolution of the SFR 
maps is significantly improved. Despite the presence of a 
few artifacts - caused e.g. by the detector 16 gaps, or by the 
imperfect superposition between tiles - a few features are 
remarkable in the maps of Fig. |13| 


• populations of all ages older than log(t/yr) >9.7 present 
a nearly-elliptical distribution of their stars; 

• the elliptical shapes persist for ages between 9.1 < 
log(t/yr) <9.7, although at these ages the distribution seems 
to have a sharper decline in density at the NE end of the 
ellipse; 

• at younger ages, however, even more significant distor¬ 
tions appear, especially at the SW side of the SMC, which 
starts to clearly delineate the Wing for all ages log(t/yr) ^ 
8.3; 

• finally, at very young ages (log(t/yr) < 8.3), a very 
patchy distribution of star formation appears, more concen¬ 
trated in the SMC bar and Wing. A prominent and wide 
episode of star formation is evident at the NE extremity 


of the bar at the log(t/yr) = 7.4 age bin. This region is 
well known from previous optical surveys and contains the 
prominent young cluster NGC 346. 


The presence of young stars in the Wing starting at 
log(t/yr) < 8.3 (t < 0.2 Gyr), is consistent with previous 
simulations which are based on the LMC-SMC-Galaxy in¬ 
teraction ( Yoshizawa V Noguchi|2003| Yozin V Bekki|20T4 ). 
These simulations clearly show that the formation of young 
stars in the Wing is due to the last LMC-SMC interaction 
about 0.2 Gyr ago, when the outer gas disk was tidally com¬ 
pressed during the very strong tidal LMC-SMC interaction 
(with the distance of the two being ~ 10 kpc). Our obser¬ 
vations appear to confirm the LMC-SMC-interaction origin 
of the Wing formation. 

Our maps can be compared with the large-scale popula¬ 
tion maps derived by Harris V Zaritsky (2004) from MCPS 
data. Such maps were obtained with a very different resolu¬ 
tion, and probably reflect many of the differences in the data 
quality and implementation of the SEH-recovery method. 
Although the comparison is not easy, some clear similarities 
are evident at first sight, especially in regard to the distri¬ 
bution of the young star formation. 

However, an aspect worth mentioning is that we do not 
find evidence of the large ring of star formation found by 


Harris V Zaritsky (2004) at ages of 2.5 Gyr, suggesting that 


this feature might have been an artifact derived from their 
significantly shallower photometry than ours, added to the 
uncertainties intrinsic to the optical work - such as for exam¬ 
ple the larger sensitivity to differential extinction of MCPS. 

These maps compare well with those inferred from the 
distribution of classical Cepheids (tracing the young, t < 
100 Myr, population) and RR Lyrae stars (tracing the old. 


10 Gyr, population) in the SMC. This is shown in Eig. 14 


where cyan points are classical Cepheids and black points are 
RR Lyrae stars. Although OGLE HI and EROS-2 cover only 
part of the area analyzed in this paper (specifically, entirely 
tiles: SMC 5_4, 4_4, 5_3, 4_3 and 3_3, half of the tile 4_5, and 
~ 1/4 of tiles 3_5 and 6_3) it is quite clear that the density 
of classical Cepheids peaks in tiles SMC 5_4, 4_4 and 4_3, 
which are dominated by the young population, according to 
our SER(t) maps. On the other hand, the RR Lyrae stars 
are distributed rather smoothly all over the SMC area - 
although in more central SMC tiles many of them probably 
do not show up because “buried” amid the overwhelming 
young population. 


© 2013 RAS, MNRAS 000, [Tpi] 
































16 Rubele et al. 


SMC' ttKiss 6.yono 
-70 n-. . .— 1 -n- 



28.5 27 25.5 24 22.5 21 10.5 1« 16.5 15 15.5 12 li:i.5 9 

RA Idegl 



28.5 27 25.5 24 22.5 21 19.5 18 16.5 15 15.5 12 111.5 9 

RA Idegl 


SMc; tiiHxs 7.8onn 



1200 

1000 

800 

600 

400 

200 

0 


SMC'tmss 8.1000 
-70 n-1—'-■ I ■ I -1—---I-1—' 



28.5 27 25.5 24 22.5 21 19.5 IS 16.5 15 15.5 12 111.5 9 

RA Idegl 


SMC' ttKiss 8.5000 
-70 n-1—■-■ I ■ I-1—---I-1—' 



2^5 27 2,5.5 24 22.5 21 19.5 IS 16.5 15 15.5 12 111.5 9 

RA Idegl 


SMC' tiiHxs 8.,'iOOO 



SMC' ttKiss 8.7000 
-70 n-1—'-■ I ■ I-1—---I-1—' 



28.5 27 2,5.5 24 22.5 21 19.5 IS 16.5 15 1,15 12 li:i.5 9 

RA Idegl 



-70 

70.4 

-70.8 

-71.2 

-71.6 

'B.jia 

u 

2J72.8 

G 

73.2 

-73.6 


28.5 27 2,5.5 24 22.5 21 19.5 IS 16.5 15 15.5 12 111.5 9 

RA Idegl 



28.5 27 25.5 24 22.5 21 19.5 IS 16.5 15 15.5 12 111.5 9 

RA Idegl 


400 

200 


0 



-70 

70.4 

-70.8 

-71.2 

-71.6 

B72.4 

U 

4J72.8 

G 

73.2 

-73.6 


28.5 27 2,5.5 24 22.5 21 19.5 IS 16.5 15 15.5 12 111.5 9 

RA Idegl 



-70 

70.4 

-70.8 

-71.2 

-71.6 

'B.llA 

U 

2J72.8 

G 

73.2 

-73.6 


28.5 27 25.5 24 22.5 21 19.5 IS 16.5 15 15.5 12 111.5 9 

RA Idegl 



28.5 27 25.5 24 22.5 21 19.5 IS 16.5 15 15.5 12 111.5 9 

RA Idegl 


400 


200 


Figure 15. The same as in Fig. |13| but no'w on a colour scale that reveals the total stellar mass formed in each age bin. 
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Figure |15| continued. 



a fXo (deg) 


Figure 14. Distribution of classical Cepheids (cyan points) and 
RR Lyrae stars (black points) in the SMC. There are 4630 
Cepheids from OGLE-III plus 165 from EROS-2, and 2475 RR 
Lyrae stars from OCLE III plus 133 from EROS-2. The black 
boxes are the VMC tiles analyzed in this work. 


Since the SFR(t) maps are built for bins almost-equally 
spaced in log(t/yr), they do not reflect the total mass of 
stars built in each epoch of the SMC history. To derive the 
latter, we redo the maps re-scaling them to the total stellar 
mass formed in each age bin. The results are presented in the 
maps of the stellar mass assembly of Fig.|15| As can be seen, 
these maps reveal that the young star formation constitutes 
only a minor fraction of the total stellar mass formed in 
the SMC. It is remarkable however that the log(t/yr) = 9.7 
age bin (4.0 to 6.3 Gyr, or ~ 5 Gyr) turns out to contain a 
significant fraction (34 %) of the total stellar mass formed 
in the SMC. 

This latter result is in agreement with |Rezaeikh et al.| 
(2014), who, analysing the iC-band luminosity function of 


long period variables found/classified by Soszyhski et al. 
(2011) and Ita et al.| ( |2004| , also find the main epoch of stel¬ 
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Figure 16. The SMC global SFR(t) (solid red line), together with 
its error bars. Also shown are the relative stellar mass formation 
as a function of age (green dotted line), and the cumulative stellar 
mass normalised to the total formed mass (of 3.85 x 10® Mq; red 
dotted line). 


lar mass assembly in the SMC to have occurred at around 
6 Gyr ago, at a rate of ~0.28 M© yr~^. 


4.5 The global SFR(t) 

With the detailed SFR(t) and population maps to hand, we 
can perform a rough comparison with the SFR(t) derived by 
other authors. Before that, we build the SMC global SFR(t) 
as illustrated in Fig. |16| Considering the large coverage of 
the central parts of the SMC provided by the present VMC 
data, it is likely that this global SFR(t) describes most of 
the stellar mass formed by this galaxy. 

Figure shows a few periods of more intense star for¬ 
mation in the SMC. Essentially, it tells us that the SFR(t) 
was relatively modest at very old ages, with mean values be¬ 
low O.O15M0 yr“^. SFR intensified significantly at an age of 
~ 5 Gyr, reaching a peak of O.ObM© yr“^; during this period, 
a large fraction of the stellar mass of the SMC was formed. 
At intermediate ages, the SFR(t) declined to values as low 
as O.O25M0yr“^, but with another peak of O.O4M0yr“^ 
at ~ 1.5 Gyr. Finally, more age resolution is provided at 
young ages, for which it is apparent that the very recent 
star formation has increased by factors between 2 and 5 
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with respect to the values dominant at intermediate-ages. 
The young SFR(t) peak is found in the second age bin, span¬ 
ning from 14 to 40 Myr. Figurealso shows the cumulative 
stellar mass formed since the oldest ages. It turns out that 
3.85 X 10^ Mq of stars were formed in the area presently 
sampled by VMC; this estimate obviously depends on the 
adopted initial mass function, which comes from |Chabrier| 
( 20 ^. 


Comparing our results with the literature is not easy, 
considering the significant number of studies in the past, 
based on very different datasets and methods. The era of 
quantitative SFH derivation starts in practice with a small 


region with available HST/WFPC2 data analysed by Dol¬ 


phin et al. (2001). It was followed by the impressive work 
by Harris & Zaritsky ( |2004| , based on the MCPS, which 
marked the field for years. Although the MCPS resolution 
was not comparable to the HST one, it had the great advan¬ 
tage of completely covering the SMC main body. The main 
findings from [Harris fc Zaritsl^ (2004) were the presence of 
several periods of enhanced star formation, with one in par¬ 
ticular (at 2.5 Gyr) causing a large-scale annular structure 
in the SMC. As already mentioned, we do not find evidence 
of such a structure. 

In more recent years, different sets of deep photomet¬ 
ric data were used for the same goal of deriving the field 
SFH (e.g. Chiosi Vallenari|2007| Noel et al.|2009 Cignoni 


et al.|201^ Weisz et al.|2013 ); although they have reached 

a general consensus regarding some basic SFH features of 
the SMC, they have also been dissimilar to some extent. For 
ages > 12 Gyr, the SFH seems to be consistent with a rela¬ 
tively low star formation, as judged from the CMD analyses 
of deep HST photometry of seven spatially diverse fields by 


Weisz et al. (2013), as well as with ground-based BR pho¬ 


tometry of twelve distinct fields spread across the galaxy by 


Noel et al. (2009). According to Noel et al. (2009), such a 


low old SFR is similar at all radii and azimuths. No region 
inside 4.5 kpc from the SMC centre has been found to be 
dominated by the old population. These findings by [Noelj 
et al. (2009) and Weisz et al. (2013) broadly agree with our 


conclusions, that the SFR has been modest at all ages larger 
than 5 Gyr. 

Relatively old and intermediate-age peaks have also 
been found in the SMC SFH at ages of ~4.5 and 9 Gyr, and 
suggested to be the result of tidal interaction with the LMC 


(Weisz et al. 2013). In the western fields studied by Noel 


et al. (2009), the older enhancement splits into two peaks, 
located at ~ 8 and M2 Gyr. However, HST VI CMDs of 
six other fields located at projected distances of 0.5 to 2 de¬ 
grees from the SMC centre resulted to exhibit a slow star 
formation pace since the galaxy formation until ~ 5 Gyr 
ago, when the star formation activity started to increase 
sharply (Cignoni et al. 2013b). Such an increase peaks at 


~ 1.5 — 2.5 Gyr. Additionally, the level of the intermediate- 
age SFR enhancement systematically increases towards the 
centre. Whether gas from the outer regions has been cen¬ 
trally tunneled still remains an issue, due to uncertainties 
in the geometry of the SMC ( Weisz et al.||2013 ). As far as 
we are aware, the conundrum about the existence of sev¬ 
eral peaks at ages older than 2 — 3 Gyr depended on the 
unavailability of datasets for more areas across the SMC. 
In regard to these issues, our data clearly confirms (1) the 
peaks of star formation at ages of 1.5 and 5 Gyr, and (2) 


that the peak of SFR at 1.5 Gyr appears, effectively, more 
concentrated than the 5 Gyr one, as can be appreciated in 
the population maps of F ig, A similar result was ob¬ 
tained by [Rezaeikh et ah ( 2014| ), who finds that the SMC 
experienced two main episodes of star formation at ~ 6 Gyr 
(log(t/yr) = 9.8) and ~0.7 Gyr (log(t/yr) = 8.8), being the 
latter more concentrated in the central parts. 

Noteworthily, the 1.5-Gyr old peak of star formation 
finds correspondence in a peak of cluster formation occur¬ 
ring at about the same time in the SMC, LMC, and the 
MW ( Piatti|2010 2011b|a ). This coincidence might be indi¬ 
cating a close encounter with the Milky Way at these ages, 
although such an event is strongly disfavoured by the most 
recent proper motion data for the LMC ( Besla et al.||20Q7 
Kallivayalil et al^|2013 ). 


As far as ages younger than 500 Myr are concerned, the 
SFH is mostly driven to the SMC innermost regions, where 
the major star formation events took place. This age range is 
featured by several peaks of star formation which have pref¬ 
erentially occurred towards the eastern side of the galaxy 
( Noel et ^[2009 ). Indu & Subramaniam (2011) also found 
evidence of a shift for the centroid of population younger 
than 500 Myr and up to 40 Myr in the direction of the LMC 
using the VI OGLE-HI and the MCPS. Our results broadly 
confirm these suggestions, but provide substantial improve¬ 
ments in the spatial resolution and area of the derived young 
SFH. 

In summary, some features of our SFH find correspon¬ 
dence with the results in the literature, others do not: 

• The initial period of low SFR is supported by the HST- 


based studies by Weisz et al. (2013) and Cignoni et al. 


2013a), and by the ground-based study by Noel et al. (2009). 


In the western fields studied by Noel et al. (2009), the 


(2012 


older enhancement is found to split into two peaks, located 
at ~8 and ~ 12 Gyr. Although our data does not extend 
as far to the west as theirs, we do not find indications of a 
similar split, since our age bins imply a coarser resolution of 


the oldest SFR(t) with respect to Noel et al. (2009). 

• The period of intense SFR at 5 Gyr finds strong support 


in the HST-based work by Weisz et al. (2013), who date this 


et al. (2012). 


event at 4.5 Gyr, and also in Noel et al. (2009) and Cignoni 


Harris & Zaritsky (2004) find bursts of star formation 


at 2 —3 Gyr, 400 Myr, and 60 Myr. The first and last of these 
periods seem to have correspondence in our data, although 
at slightly different ages. 


Cignoni et al. (2013b) find an increased SFR peaking 


at ~ 1.5 — 2.5 Gyr, which is consistent with suggestions of 
a close encounter with the Milky Way at recent time. We 
essentially find the same, at 1.5 Gyr. 


From a different kind of data, Rezaeikh et al. (2014) 


finds main episodes of star formation at ~ 6 Gyr and ~ 
0.7 Gyr, at rates in reasonably good agreement with our 
data, and with a concentration of the younger star formation 
towards the center, similar to the trends we find. 

As a final note, we point out that all the previous works 
on the SMC SFH seem complementary, in the sense that 
they analyse the SFH at different spatial scales and with 
very different sensitivities regarding the old SFH. From the 
point of view of the present work, they all suffer from the 
same problem: a higher sensitivity to the presence of differ- 
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Figure 17. The SMC global AMR (solid black line with error 
bars). The colour points show the distribution of the stellar partial 
models in the AMR diagram. The olive points with error bars 
show the distribution of clusters analysed i n |Piatti| (|2012| ). The 
violet points are the mean values obtained by |Carrera et al. f poo^ 
from the combination of 13 SMC fields; in this case the vertical 
error bars indicate the metallicity dispersion inside the age bins 
delimited by the horizontal error bars. 


ential extinction. None of these surveys can rival VMC in 
terms of spatial coverage, except for the MCPS; and yet, 
MCPS is the survey that more dramatically illustrated how 
the SFHs derived from optical data are sensible to differ¬ 
ential extinction: Harris & Zaritsky (2004) even had to re¬ 
sort to different amounts of extinction affecting red and blue 
stars, in order to improve their CMD fits. With its sensitiv¬ 
ity to dust reduced by a factor of at least three, as compared 
with a survey based on VI data, this is one aspect where 
VMC has probably made significant progress. 


4.6 The metallicity evolution 


In the context of the SFH-recovery work, the AMR can be 
treated either as an integral part of the derivation method, 
or as a by-product. The difference is in whether the AMR 
derived from field observations is used to force the SFH solu¬ 
tion to occupy a certain line in the age-metallicity space, or 
whether the best-fit SFH solution is used to indicate which 
AMR actually holds. We have used the latter approach, al¬ 
lowing the AMR to vary within a wide area in the age- 
metallicity space, providing enough freedom for slightly dif¬ 
ferent AMRs to be derived between different regions. Need¬ 
less to say, such an approach would not have been possible 
were VMC not to cover an extremely large area of the SMC. 

The several space-resolved AMRs have already been in¬ 
troduced in Fig.|^ In addition, we have produced the global 
AMR illustrated in Fig. which is simply the weighted 
mean [Fe/H] computed over ah of the SMC, for each age 
bin. The general features are clear: the [Fe/H] has generally 
increased in time, with the exception of a period of slight 
decrease in metallicity taking place for ages slightly younger 
than I Gyr. 

Several independent determinations of the AMR exist 
in the literature, most of them derived from the study of 


ages and metalhcities in star clusters. Piatti (20IIc) used 


the largest known SMC cluster sample with ages ^ I Gyr 


and metalhcities put onto a homogeneous scale, and found 
two periods of enhanced cluster formation at t ~ 2 and 5 — 6 
Gyr, taking place throughout the entire body of the galaxy. 
He also found the absense of a metallicity gradient, and a 
spread in metallicity for clusters older than ~ 7 Gyr. A com¬ 
promise between covering a relatively large area and reach¬ 
ing the oldest main sequence turn-offs was accomplished by 


Piatti (2012). He found that the held stars do not possess 


gradients in age and metallicity, and that stellar popula¬ 
tions formed since ~ 2 Gyr ago are more metal-rich than 
[Fe/H] ~ —0.8 dex and are conhned to the innermost re¬ 
gion (semi-major axis ^ 1°). He also compared the held 
star AMR to that of the star cluster population with ages 
and metalhcities in the same held scales, and found that 
clusters and star helds share similar chemical evolution his¬ 
tories. This latter conclusion seems to be fully supported 
by the comparison with our results for the held, shown in 
Fig. |I7| For comparison, the hgure also presents the mean 
AMR obtained by Carrera et al. (2008) from caltium triplet 
observations in 13 different SMC helds. Note that there is 
a good agreement between the different AMRs for all ages 
larger than 2 Gyr. For the youngest age bin, the comparison 
is meaningless because the RGB-dominated sample by |Car-| 


rera et al. (2008) does not provide age information for stars 
younger than ~ I Gyr. 

One of the most remarkable points in the global AMR 
of Fig. is the peak in metallicity at 1.5 Gyr, followed 
by a decrease at 0.8 Gyr, and then hnally followed by a 
steady increase that leads to the present-day metalhcities. 
This indicates an event of dilution in the metallicity having 
occurred ~ 1.2 Gyr ago. 

In order to better understand this feature, in Fig. |I8| we 
plot the spatial distribution of the metallicity in each sub re¬ 
gion analysed, in the relevant age interval. It is remarkable 
that the peak in metallicity manifests itself hrst in the age 
bin log(t/yr) = 9.3 (2 Gyr), in the form of a higher metahic- 
ity in the inner SMC regions. For the age bin log(t/yr) =9.1 
(1.3 Gyr), the entire SMC is observed to present a higher 
metallicity. Finally, for the age bin log(t/yr) = 8.9 (0.8 Gyr), 
the metallicity has already decreased again, but for the 
Northern outskirts of the SMC - which, coincidentally, are 
the only regions to not have formed stars recently. Although 
the dynamics of the SMC stellar populations is quite uncer¬ 
tain over such timescales, the effect appears as if a major 
phase of accretion of metal poor gas has taken place in the 
SMC history, interrupting the normal increase in the sys¬ 
tem’s metallicity, at an age of I Gyr, so that the new star 
formation took place starting with a metal-poorer gas. In¬ 
terestingly, the period of global decrease in the SMC metal- 
licity follows the I.5-Gyr peak in the global SFR(t). This 
concomitance might be indicating a major event in the SMC 
history, like for instance the merging with a gas-rich com¬ 
panion galaxy, or simply a closer passage to the LMC or the 
Milky Way. 

Indeed, this latter interpretation seems consistent with 
the simulations by Yozin & Bekki (2014), which predict that 
the LMC-SMC-Galaxy interaction could have triggered en¬ 
hanced star formation in the SMC about 1.5 Gyr ago. This 
is also the epoch of the Magellanic Stream formation. They 
hnd that the I.5-Gyr ago interaction can strongly distort the 
entire SMC body so that a large amount of gas in the outer 
part (thus metal-poor) can be funneled to the central bar 
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Figure 18. Spatial distribution of the mean metallicity derived for each age bin between log(t/yr) 


8.9 and 9.3 (0.8 to 2 Gyr). 


region, where enhanced star formation can occui[^ When 
the outer metal-poor gas is transferred to the central region 
to be consumed by star formation, the new stars should 
have lower [Fe/H] initially. However, owing to the chemical 
evolution after the burst, [Fe/H] should increase again. A 
subsequent study by Yozin & Bekki (2014 see their figure 
14) does not show clearly the dip in the AMR, perhaps be¬ 
cause they averaged the AMR for the entire region. In the 
present work, however, the same galaxy-wide average still 
causes the metallicity dip to be clearly seen in the data. 

On the other hand, as a cautionary note, we recall 
that the > 1.5 Gyr interval also corresponds to the ages 
at which the RGB quickly develops in stellar populations, 
so we cannot exclude that our method reacts to this sharp 
evolutionary transition by introducing some artificial fea¬ 
ture in the AMR at those ages. However, in support of the 
“galaxy-interaction” hypothesis, is that most of the massive 
intermediate-age clusters in the LMG and SMG turn out to 


be observed exactly at ages close to 1.5 Gyr (see e.g. Gi- 


rardi et al. 2013). Moreover, the most massive LMG and 


SMG clusters of ages between 1 and 2 Gyr are known to 


present multiple turn-offs and dual red clumps ( Goudfrooij 
et al.|200^ Girardi, Rubele Kerber|2Q09 ), which might be 


related to the same phenomenon: i.e., a galaxy interaction 
that has disturbed the ISM enough to drive the accumula¬ 
tion of gas in the cluster interiors, and the subsequent star 
formation there. Although this might sound as rather spec¬ 
ulative, this picture could explain a wide range of observa¬ 
tional facts, without resorting to extreme assumptions. 


5 SUMMARY AND CLOSING COMMENTS 

In this work, we have provided SFH, distance and extinction 
maps across an area of 14 deg^ of the SMC. The maps are 
computed with an initial resolution of 20^ which is the size of 
the 120 nearly-square subregions used in the SFH-recovery. 
Although maps of similar or better resolution have been 
presented before, our work is distinguished from them by 


^ Such simulations do not predict enhanced star formation at 
~ 5 Gyr, because this is beyond the scope of the currently exist¬ 
ing LMG-SMG interaction models. A possible burst (via merger) 
around 7.5 Gyr ago in the SMG is discussed by |Tsujimoto fc Bekki| 
( |2009| >. The 5 Gyr burst could constitute a useful constraint on 
such models, being described e.g. as the epoch of the first arrival 
close to the MW, or the epoch when the SMG started its strong 
interaction with the LMG. 


basically three crucial aspects: (1) the use of uniform and 
high-quality near-infrared photometry, reaching stars near 
the oldest main sequence turnoff across all of the observed 
area and with a reduced sensitivity to differential extinction; 
(2) the simultaneous derivation of all the quantities from the 
best global ht of the GMDs; (3) the extremely wide area that 
has been analysed, in a homogeneous way. 


These results already allow us to revise previous results 
regarding the SMG morphology. We hnd that it is system¬ 
atically tilted in the E-W direction, with a mean inclination 
angle of 39°. However, there are enough signihcant displace¬ 
ments with respected to an inclined disk, to allow us to 
infer that the SMG is signihcantly distorted and similar to a 
warped disk. In addition, a large distance spread along the 
line of sight is clearly detected in the SE end of the SMG, 
which corresponds to the end of its Wing. After deriving the 
SEH for ah these subregions, we assign to every observed star 
a probability of belonging to a given age-metallicity interval. 
This allows us to make high-resolution maps of the several 
stellar populations observed across the SMG. We can, for in¬ 
stance, illustrate the flocculant nature of the very young star 
formation, and identify the nearly smooth features drawn by 
the older stellar generations. We identify the centre of the 
intermediate-age to old stellar populations to be located at 
coordinates (ac = 12.60°,(5c = —73.09°), which is 3.7° dis¬ 
placed with respect to the kinematical SMG center deter¬ 
mined by Stanimirovic, Staveley-Smith & Jones (2004). We 
document the formation of the SMG Wing at ages younger 
than ~ 0.2 Gyr, and the dramatic features drawn by the 
young star formation in the SMG Bar, which peaked at ages 
of about 40 Myr. Even more importantly in terms of the 
SMG history, we clearly detect periods of enhanced star for¬ 
mation at ages of 1.5 and 5 Gyr ago. The former is probably 
related to an interesting feature we hnd in the AMR, that 
points to the ingestion of metal-poor gas taking place at ages 
slightly larger than 1 Gyr. The latter, instead, constitutes a 
major period of stellar mass formation in the SMG. Indeed, 
we conhrm that the SER was moderately low at even older 
ages in the SMG. 


These hndings are still to be further explored, and fully 
interpreted in terms of the possible scenarios for the evolu¬ 
tion of the SMG. Although they clearly conhrm many pre¬ 
vious claims in the literature, they also provide signihcant 
more detail on several aspects, such as the wide-scale distri¬ 
butions of the populations with different stellar ages. 

The SEH data are available at http://stev.oapd. 
inaf.it/VMC_SFHdata/ The analysis of additional VMG 
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data, that over the next four years will complete the SMC 
and LMC mosaic, is in progress. 
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APPENDIX A: HOMOGENIZED DEEP TILES 
AND PSF PHOTOMETRY 


To generate a deep tile we combine images from different 
epochs (at least three in Y and J and twelve in the 
band), each one formed by the combination of six pawprint 
images. These images are stacked single exposures reaching 
a total exposure time of 400 seconds in Y and J, and 375 
seconds in iFg, for each pawprinlj^ Variations of the seeing 
occurring over these timescales could affect the PSF shape 
on the final deep tile, as a function of position. To avoid this, 
we apply a homogenization method to correct the PSF on 
each detector, in each pawprint image, in each epoch, to a 
constant and homogeneous PSF model, before generating a 
deep tile. We proceed as follow: 


(i) First we derive the effective PSF model (EPSF) on 
pawprint images of a given epoch. To do that we assume 
a constant PSF across each detector of the VIRCAM cam¬ 
era, and use the IRAF/DAOPHOT tasks to find the best-fit 
PSF model. Usually, it corresponds to a Moffat function pro¬ 
file. Examples of such EPSFs are shown in the left panel of 
Fig.|AT| 

(ii) Afterwards, we generate a symmetric Moffat function 
reference PSF model (RPSF) with a half-flux radius (HFR) 
equal to the largest HFR of the EPSF model of all detectors, 
in all pawprints images, at a given epoch. 

(iii) Subsequently, using a Fourier deconvolution method 
we derive the kernels that, when convolved with the paw¬ 
print images, convert them to images with a homogeneous 
RPSF model. An example of kernels and its RPSF model, 
applied on a pawprint image are illustrated in the central 
and right panels of Fig. 1^ 

(iv) Finally we combines all homogenized pawprint im¬ 
ages with the SWARP tool ( Bertin et al.|200'2 ) in the same 
way as described in Rubele et al. (2012), thus generating a 
deep tile image with a homogeneous PSF. 


With this method, we obtain epoch images (in both 
pawprint and tile formats) with a homogeneous PSF, of 
course degraded to the worst pawprint at a given epoch. 


® The average observing time across the deep tile is about twice 
these values. 



k*vsa 



Figure A2. Difference between the PSF photometry in T, J, 
and Ks (panels from bottom to top, respectively) and the VSA 
aperture photometry, for the tile SMC 3_5, as a function of mag¬ 
nitude. 


but gaining on quality of photometry in the hnal deep tile. 
The hnal deep tile image in fact presents a constant, round 
and homogeneous PSF, that allows us to perform aperture 
or PSF photometry on the entire tile without correction on 
the hnal catalog and optimising to improve the PSF pho¬ 
tometry. Usually, variations in the seeing less than 35%, that 
represent most of seeing differences in the VMC data set, are 
well corrected with this methodology and the degradation 
of the hnal deep tile is less than 20%, compared with the 
original (non-homogenized) deep tile. 

Having a homogenized tile image is also important for 
the evaluation of errors and completeness through very ex¬ 
tensive artihcial star tests (AST, see Sect. 2.2), preventing 
the occurrence of systematic errors due to differences be¬ 
tween the PSF model used to simulate the artihcial stars, 
and the PSF of real stars on the image. 

Finally, a comparison between our homogenized PSF 
photometry and the VSA aperture photometry is shown in 
Fig. |A2| It reveals that the two photometries are similar, 
to within 0.1 mag, down to depths of 22 mag in Y and J, 
and 21.5 mag in K^. Of course, differences are larger in high 
density helds or in the centres of star clusters, where the 
PSF photometry is expected to be more reliable. 

Figure |A3| shows an example of the efficiency of the 
DAOPHOT parameter “sharp” in separating between point¬ 
like sources and extended objects - supposed to be stars and 
galaxies, respectively. Simple cuts in the sharp versus mag¬ 
nitude diagram (central and right panels) are able to isolate 
most stars, which have sharp values close to zero. Significant 
groups of objects with different sharpness are found both 
at the brightest and faintest magnitudes (black and cyan 
points). The former are simply partially-saturated stars at 
Ks < 11.4 mag. The latter are mainly background galaxies. 
Examination of the CMDs (left panel) reveals that they are 
mainly located in the red part of the CMD. Fortunately, the 
bulk of the galaxies are redder than most SMC and MW 
foreground stars. 
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Figure Al. Left Panel: An example of the derived PSF model on each detector on a pawprint image at a given epoch. Central panel: 
The deconvolved kernels for the same detectors. Right panel: The RPSF model used to produce the kernels. 



Figure A3. Left panel: CMD of tile SMC 3_5 showing the position of point-like objects (violet), and extended objects (black and 
cyan) selected using the sharp parameter in Kq and Y bands, respectively. Central panel: the sharpness parameter as a function of the 
iFg-band magnitude, with yellow dots marking objects selected as point-like and black dots marking the extended ones. Right panel: 
the same for the Y band, with burgundy and cyan dots for point-like and extended objects, respectively. 


APPENDIX B: PHOTOMETRIC ZEROPOINTS 


As described in Rubele et al. (2012), we have estimated the 


differences between the calibration equations provided by 
CASU (eq. |B1| and those predicted by the theoretical mod¬ 
els which are, by construction, strictly on a Vegamag system 
(where Vega star has a null magnitude in all filters). The 
same process is repeated for this work, since we are using 
a different version of VS A data, and improved photometry. 
The basic idea is to reproduce, with the aid of stellar mod¬ 
els, the same sort of photometric data used by CASU in 
the original calibration of VISTA data from 2MASS obser¬ 
vations. 

We proceed in the following way: First we generate a 
synthetic model of the Milky Way (MW), containing both 
2MASS and VIST A filters, using the TRILEGAL code ( [gT 
rardi et M]|2QQ5 ). The simulated region has a total area of 


1.0 deg^ towards the SMG. The typical 2MASS photomet¬ 
ric errors were added to the simulations following [Bonatto^ 


Bica V Girardi (2004). They dominate the error distribution 


compared to VISTA data. Then, we fit the distributions of 
these points in the colour-colour diagrams shown in Fig. |Bl| 
keeping the slope fixed at the same value as adopted by 
GASU, and adopting a sigma clipping of 3cr. For the best-fit 
equations we obtain: 


Avista — e/2MAss 
VviSTA — e/ 2 MASS 
JviSTA — J 2 MASS 
i^VISTA — i^2MASS 
VISTA — ^s2MASS 


1.025(J-if)2MAss + 0.082 
0.610(J-if) 2 MAss + 0.074 (Bl) 
-0.077(J-i7)2MAss + 0.026 
0.032(J-i7)2MAss + -0.017 
0.010( J—iFs)2MAss + 0.003 


These fits imply offsets of 0.082 mag in Z, 0.074 in 
y, 0.026 in J, —0.017 in H and 0.003 in Ag, between the 
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Figure Bl. Calibration to a Vegamag system. The points show 
the difference between VISTA and 2MASS magnitudes versus the 
2MASS colours of simulated Milky Way model stars on an area 
of ~ 1 deg^. The simulated stars are generated assuming the 
observed distrbution of 2MASS errors. Green points show the 
distribution of stars inside the adopted 3cr clipping limit (black 
dotted lines). The black dashed dotted line shows the best-fit 
linear relation, with slopes fixed as in Eqs. |B1| 


converting them to the same set of ZPs as the observations, 
before starting the work of SFH- recovery]^ 
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